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We construct a minimal model that naturally generates small neutrino masses and provides a 
dark matter candidate. The symmetry making the dark matter particle stable simultaneously sup- 
presses neutrino masses to appear first at 3-loop level. This provides an elegant explanation for 
the observed mass hierarchy ruv/v ~ 10"^"^, without introducing right-handed neutrinos. The dark 
matter particle is part of an inert scalar doublet that plays a crucial role in the radiative neutrino 
mass generation. The model gives distinct predictions: normal neutrino mass hierarchy, a flavour 
mixing angle ^13 ~ 7° — 11°, new charged scalars at around the electroweak scale and a dark matter 
candidate in the mass range 50 to 75 GeV. 
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I. INTRODUCTION. 

The existence of non-baryonic dark matter in the Uni- 
verse and the experimental observation of nonzero neu- 
trino masses may be regarded as the most direct and 
compelling evidence of particle physics beyond the stan- 
dard model (SM). However, both the origin of neutrino 
masses and the nature of dark matter are still unknown. 
A scenario being able to incorporate both phenomena in 
a unified framework would then be very attractive. 

One of the best-motivated dark matter scenarios is that 
of stable weakly interacting massive particles (WIMPs), 
produced as a thermal relic from the early Universe. 
Among the simplest realizations of this WIMP scenario 
is the inert doublet model [l|, . The SM is extended by 
a scalar doublet $2, and the dark matter scalar is made 
stable due to an exact Z2 symmetry under which the new 
field has odd parity. The inert doublet model is currently 
constrained by dark matter searches as well as by parti- 
cle collider data, but still a large region of its parameter 
space is allowed. 

From the perspective of neutrino physics, currently the 
most popular way to generate small neutrino masses is 
the see-saw mechanism (see 0] for a review). In its sim- 
plest variant, it postulates the existence of very massive 
S'C/(2)l X U[1)y singlet right-handed neutrinos. Other 
realizations involve the existence of very heavy scalar 
triplets or fermionic triplets. Although elegant, this 
mechanism is difficult (if not impossible) to test experi- 
mentally, as the masses Mi of the new states are typically 
much larger than the electroweak (EW) scale. Mi ^ v. 

Small neutrino masses can also be generated via ra- 
diative corrections. This has been explored with explicit 
lepton number violation in extensions of the scalar sec- 
tor of the SM. As opposed to the see-saw mechanism, this 
approach generates small neutrino masses without rely- 
ing on a very high energy scale. One simple realization 
of this idea is the Zee model [1], in which the SM field 
content is enlarged by a second scalar doublet $2 and 
a charged scalar singlet 5'+. Another simple scenario is 



the Babu model which replaces the scalar doublet 
<I>2 in the Zee model by a doubly charged singlet scalar 
field p"*"*". Other realizations of this idea, with or without 
right-handed neutrinos, have also been proposed [l|,@-[ll- 

However, these scenarios of radiative neutrino mass 
generation in Q and @ (together with many others, such 
as 0, [1] ) cannot at the same time contain a viable dark 
matter candidate. A stabilizing symmetry for some of the 
new fields would in fact forbid the very terms responsible 
for the generation of neutrino masses (see however [l|, 0| 
for interesting scenarios including right-handed neutri- 
nos and dark matter particle candidates). This is for 
example the case in the Zee model if an odd Z2 parity 
is assigned to $2- Finding a common scenario for radia- 
tive neutrino mass generation and dark matter is then a 
nontrivial task. 

In this letter we construct a minimal model, that gen- 
erates neutrino masses radiatively and provides a stable 
dark matter candidate, via an extended scalar sector with 
an exact Z2 symmetry. We do this by adding to the SM 
two scalar singlets, p'^^ , , and a scalar doublet $2- 
The fields 5+ and $2 have odd Z2 parity (while other 
fields do not transform under this symmetry), and there- 
fore a variation of the mentioned inert doublet model 
of dark matter is automatically embedded into the sce- 
nario. Due to the Z2 symmetry and the field content 
of the model, Majorana neutrino masses are first gener- 
ated at the 3-loop level, naturally explaining the large 
hierarchy rriy/v ^ 10"^'^ as due to the loop suppression 
(g2/l67r2)3 ^ 10-13 [g being a EW-sized coupling). This 
scenario then provides an intrinsic and interesting link 
between the stability of the dark matter candidate and 
the smallness of the neutrino mass scale. 

The letter is organized as follows: in section II we in- 
troduce the model, obtain the neutrino mass matrix and 
verify that it produces the observed pattern of neutrino 
masses and mixings. In section HI we discuss the re- 
quirements on the dark matter candidate. In section IV 
wc analyze the constraints on the model coming from 
lepton flavour violating (LEV) processes, EW precision 
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FIG. 1: The Cocktail Diagram 



tests (EWPT) and collider searches, and we comment on 
possible consequences for neutrinoless double beta de- 
cay (Qvfip). We then briefly discuss future detection 
prospects, before concluding. 



II. A MODEL FOR NEUTRINO MASSES. 

In addition to the SM fields, the model includes two 
S'C/(2)l singlet scalars (singly and doubly charged) 5"+ 
and and a scalar doublet $2- We introduce a Z2 

symmetry under which the <f>2 and 5"*" fields are odd, 
whereas p^'^ and the SM fields are even. The Z2 sym- 
metry should be unbroken after EW symmetry breaking, 
so that the lightest Z2-odd state remains stable and can 
provide a dark matter particle candidate. Given the sym- 
metry and particle content of the model, the lagrangian 
will include the following relevant terms leading to lepton 
number violation 



A/: 



Kl$2«cr2$lS' +K2P++5' S' 



As 



(1) 



The SM scalar doublet $1 and the inert scalar doublet 
$2 can in the unitary gauge be written as 



$1 = 



1 







J_f A+ 



(2) 



where v ~ 174 GeV is the vacuum expectation value of 
$1. After EW symmetry breaking, and for ki 7^ 0, the 
charged states and 5*+ will mix (the mixing angle 
being /3), giving rise to two charged mass eigcnstates 



3/3 



-C/3 



A+, 



cpS+-s,3A+, (3) 



with sp,cp = sin/3,cos/3 respectively. 

The lagrangian in Eq. ([T|) breaks lepton number explic- 
itly by two units Q, which generates a Majorana mass 



for the left-handed neutrinos. The Z2 symmetry pre- 
cisely forbids all terms that would have generated neu- 
trino masses at either 1 or 2-loop order, and therefore 
the leading contributions to neutrino masses appear at 3- 
loops through the ^Cocktail Diagram^ shown in Figure [T] 
In the basis where the charged current interactions are 
flavour-diagonal, the charged leptons e,/i,T being then 
mass eigenstates, and after summing up the contributions 
from the six different finite 3-loop diagrams in Figure [T] 
(coming from H^2j ^0 and Hq running in the loop), the 
Majorana neutrino mass matrix reads: 



™a6 — ^ab Xa Xb S^jj 



I. 



A 



(16^2)3 

where S2/3 = sin(2/3), Xa = rria/v for a = e, /i, r, and 



A = 



Pa M4 



(4) 



(5) 



The factor I^, is a dimensionless 0(1) number emerging 
from the 3-loop integral after all generic factors have been 
factorized out. Its exact value depends on the specific 
mass spectrum, and we have estimated its value using 
the numerical code SecDec flfl. The reduced masses are 



no 



and 



The dependence of mj^j, on the mass differences Atoq 
?r and Ami 



ifj^ — rrifj^ signals a GIM-like 
mechanism at play in Eq. which can be easily under- 
stood noticing that AmQ oc A5 and Am^ oc ki. In the 
limit A5 ^ the lagrangian in Eq. ([1]) conserves lepton 
number and no Majorana neutrino mass can be gener- 
ated, while in the limit ki 0, the leading contribution 
to mj^jj will appear at a higher loop order. 

We now analyze the ability of the model to reproduce 
the observed pattern of neutrino masses and mixings. 
The standard parametrization for the neutrino mass ma- 
trix in terms of three masses mi^2,3j three mixing angles 
di2, 023 J 013 and three phases S, ai, a2 reads 



m 



U U with TO 



D 



Diag(TOl,TO2,TO3) (6) 



C/ = Diag (e'"l/^e'"^/^l) 



C13C12 -C23Sl2-S23Cl2Sl3e 
C13S12 C23Cl2-S23Sl2Sl3e'* 

iSise""^ S23C13 



S23Sl2-C23Cl2Sl3e 
-S23Cl2-C23Sl2Sl3e 
C23C13 



iS 



with 



, ) and 



cos{9ij). A global fit to 



neutrino oscillation data after the recent measurement 
of 01 3 (see for example [ll|) gives Arn^i = TO2 — = 



7.62 

10- 

^23 



713 
+0.19 
0.19 



X lO-^'eV^ 



^eV" 



43+0 03 

^■^'^-0.03 



"311 = 
hO.016 
-0.017' 



'13 



71^1 — 2.55j]Q_Qg X 
0.025l°:™i^, and 



0.320 

(0.6lt^;J]^) if in the first (second) oc- 



tant for 623- Neutrino oscillations are not sensitive to 
the Majorana phases ai and a2 nor to the absolute neu- 
trino mass scale, while the value of the CP phase S is 
beyond current experimental sensitivity. In the inverted 



hierarchy scenario {Arn^^ < 0) these experimental data 
lead to |mgg| > 10^^ eV, which cannot be accommodated 
by Eq. (jl]) due to the ^ 10~^ suppression. Thus, the 
model predicts a normal hierarchy pattern for neutrino 
masses. 



The fact that the entries : 



in Eq. (j4]) are paramet- 
rically much smaller than the rest (being proportional to 
and XeX^) results in an approximate neutrino mass 
texture of the form rn^g ~ 0, ^ 0. For given val- 



ues of Af 



'31' ''12 



and 



'23' 



the four constraints 



Re [m^g] , Im [m^g] ~ and Re ["ig^] , Im [m-g^] ~ can 
in fact only be satisfied over a certain range of [1]. 



If Ar 



'21 ' 



'31' ^12 



and 



are taken as the central 



values of the global fit in then the model predicts 
0.017 < s?3 < 0.033 (0.010 < sfg < 0.012) for 023 < 7r/4 
(6*23 > 7r/4). It thus gives a related prediction for 6*13 
and the deviation of ^23 from 7r/4 (maximal mixing) to- 
wards lower or larger values. Moreover, within the al- 
lowed range for 5^3, the above constraints lead to a spe- 
cific prediction for mi, ai, a2 and 5 (as an example, 
for central values of Am|j^, Arngj^, s^jj S23J with 

^23 < 7r/4, this corresponds to mi ~ 0.0053 eV, ai ~ 88°, 
a2 ~ -70°, 5 ~ -100°). 

For the mass texture discussed above, the experimen- 
tal pattern of neutrino masses and mixings results in a 
neutrino mass matrix with the structure |mggg^| ~ 0, 
|m-J :^ 8 - 10 X lO'^ eV, | m'^^^^,^^, \^2-ix lO'^ eV. 
Compared to Eq. @ we find that in order to radiatively 
generate mjf^ and mj^^ entries of the right size (bounds 
on mj^^ and m'^^ are much weaker) wc need 



GeT4pA> 1.5 TeV 



Ca'm s^^^>0.3TeV 



(7) 



III. DARK MATTER. 



When the lightest Z2-odd states is electrically neutral 
the model has a WIMP dark matter candidate. For the 
remaining of the paper this particle will be assumed to 
be (taking ^0 would be equivalent). This WIMP sce- 
nario resembles the inert doublet model [3 and should 
share much of its phenomenology (see e.g. [l2| and refer- 
ences therein). 

The relic abundance of Hq is determined by its annihi- 
lation rate at freeze-out. In the mass range 50 - 75 GeV 
or above 520 GeV [l3| for m^fj,, the correct dark 
matter abundance can be achieved while being compat- 
ible with existing bounds from LEP, EWPT, and dark 
matter direct and indirect searches. The lower WIMP 
mass range allows to simultaneously generate neutrino 
masses of the right size in the model. The correct dark 
matter abundance may be reached by 

(a) annihilations into fermions via resonant SM scalars 
when rriHa ™h/2. 

(b) co-annihilation with cither Aq or if the mass 
splitting to Hq is of a some GeV. For fairly large 
mass splittings Am^ and Amp, co- annihilations Hq- 




FIG. 2: 1(7 (Red/Dark) and 2cr (Green/Light) allowed regions 
for AT as a function of ruHi and m/^j, for P = 7r/4, ruAo = 
250 GeV, tuho = 70 GeV and = 1 TeV. 



Aq is strongly suppressed, while co-annihilations Hq- 
Hi may still be possible, 
(c) the WIMP mass approaching my/^ where the close- 
ness to the WW threshold regulates the annihilation 
rate at freeze-out. 



IV. EXPERIMENTAL CONSTRAINTS. 

Direct searches for the doubly charged scalar /9++ at 
LEP via single production and decay into same-sign di- 
Icpton final states set a lower bound rup > 160 GeV [lH| 
(which however depends on the value of Cee)- Bounds 
from virtual p"*""*" exchange in Bhabha scattering lead to 
Gee ^ 9-7 X 10-6 GeV"^ m^ [ll,[3. More stringent lim- 
its from direct searches at Tevatron and LHC are more 
subtle to derive (as opposed to doubly-charged states 
A++ coming from scalar triplets, p"'""'" does not couple to 
W bosons). ATLAS searches for via pair production 
at LHC set a limit nip > 400 GeV [13 ■ For the charged 
states H^2: the LEP bound on chargino searches can be 
translated into an approximate bound mjj+ > 70 — 90 
GeV (depending on rriHg) . Moreover, LEP excludes 
models with m^^ < 100 GeV if iriHa ^ 80 GeV and 
mAo -mHo > 10 GeV [H]. 

The new inert fields also contribute to EWPT observ- 
ables, such as the oblique parameters S, T and U . For 
rrih — 126 GeV, the most important constraint is given 
by AT e [-0.04,0.12] at 95% C.L. [U (contributions 
to S and U arc found to be negligible). The one-loop 
contribution to T from the new fields is given by 



AT = 



1 



cfs {FHt,Ha + Fhi.Ao) 



~sl {Fh.2.Ho + Fh-2.Ao) - 2cIsIFhi,h2 - Fho.Ao 



(8) 



where K 



2 



■In- 



and is the Wein- 



berg angle. EWPT constraints can be satisfied for a wide 



4 



range of masses and mixing angles j3, and as opposed to 
the inert doublet model the present scenario allows 
for large mass splittings (see Figure [2]). 

In addition, the doubly charged scalar p++ mediates 
LFV tree-level processes such as /i -^> 3e and r — > 3e, 3//, 
and 1-loop processes like /i — > 67 and t — > 57, /i7. These 
constrain the allowed values of Cgb as a function of to^, 



the most stringent bounds being 



/i — 3e 



r - 
T - 



3e 
3/i 

-> /i+ 

e+ e~ /i~ 
e+ /i^ /.r 
67 



CeM C'ee < 2.3 X 10-"GcV"^ TO^ 
CerCee < 1.0 X lO-^GcV- ...p 



m 



C^T Cee < 8.0 X 10-^GeV"^ 
Cer Ce^ < 6.0 X IQ-^GcV"^ 
Cer C^^ < 8.0 X lO-^GcV^^ ml 
E ^^ip ^^ie < 4.7 X 10-9GeV-2 tt 



LFV constraints then favour 



> TeV, which com- 



bined with d?]) leads to large values of K2 {4,11 mp > 
K2 > mp), close-to-maximal mixing (3 ^ ■n/A and fairly 
large mass splittings Am^ , Attiq ^ . For such large 
mass splittings, satisfying the EWPT constraints then 
requires a mass spectrum > '^^Ao > ttih^ , result- 

ing in a partial cancellation of the Hi and H2 contri- 
butions to Eq. ([5]). As an prototypical example satisfy- 
ing EWPT (see Figure 2) and collider constraints, we 
take niH^ = 400 GeV, mn^ = 90 GeV, toaq = 250 



GcV, 

TeV and /3 



70 GeV and 



1 TcV, with K2 



7r/4. Neutrino masses of the right size 



are obtained for Cer ~ 0.04, C„ 



0.01, Cp 



0.001 



and Ct-t ~ 0.0001, satisfying at the same time all LFV 
bounds (together with CeeCep S 10~^ from the decay 



/i — > 3e). Branching ratios for several LFV processes 
(like — > 6+ ^~ n~ and /x — >■ ej) are predicted close to 
the current experimental bounds, and may be probed in 
the near future. 

The model may also be probed by future Ovfi/S de- 
cay experiments for not-too-small values of Cee, since 
the short-distance physics contribution to the Oz^/3/3 de- 
cay process is much larger than the one coming from m^g, 
and might actually be sizable in some cases. 

To conclude, we find it possible to minimally extend 
the SM to include neutrino mass generation and dark 
matter in a unified scenario. While giving an elegant 
explanation of the hierarchy m^/v^ the model predicts a 
nontrivial relation between ^13 and the octant of ^23, to 
be tested by future neutrino experiments. It also predicts 
large LFV, WIMP dark matter with m 50 - 75 GeV 
and new scalar states to be searched for at the LHC. 
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